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Abstract. We consider the transverse dynamical structure factor of the 
anisotropic Heisenberg spin-1/2 chain (XXZ model) in the gapped antiferromag- 
netic regime (A > 1). Specializing to the case of zero field, we use two independent 
approaches based on integrability (one valid for finite size, the other for the in- 
finite lattice) to obtain the exact two-spinon part of this correlator. We discuss 
in particular its asymmetry with respect to the 7r/2 momentum line, its overall 
anisotropy dependence, and its contribution to sum rules. 
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1. Introduction 

Since the early days of quantum mechanics, Heisenberg spin chains [I] have been 
one of the most fertile laboratories for investigating strong quantum effects in 
physically realistic condensed matter systems. Perhaps the most important conceptual 
development associated to spin chains is Hans Bethe's famous Ansatz [2] for their 
wavefunctions, leading in particular to exact thermodynamics for these and a wide 
variety of other types of integrable models (see [3J U] and references therein) . 

Another extremely important concept which studies of spin chains help to pin 
down is that of particle transmutation. Namely, when correlations become important, 
the microscopic degrees of freedom transform into effective quasi-particles barely 
resembling those of a free or weakly-coupled system. For Heisenberg spin chains, 
this occurs when spin waves (which are the exact excitations at the saturation field) 
give way to quantum solitonic excitations known as spinons [5] (which are the exact 
excitations in zero field). The power of integrability is illustrated by the fact that 
it allows this whole transmutation process to be described accurately. Since spinons 
should be understood as excitations created over the highly entangled ground state of 
the spin chain, they possess a lot of richness, at the simplest level in their dispersion 
relation, but also at a more elaborate level in their dynamics (in other words their 
ability to 'carry' correlations around in the system). The latter point is quantified by 
their contribution to various correlation functions, most important of which are the 
dynamical spin-spin correlation functions which determine inelastic neutron scattering 
amplitudes. Various guises of integrability can now be used to address the question 
of correlation functions in such strongly-correlated systems, and the purpose of this 
paper is to investigate such a case. 

Spinons come most radically to the forefront for the anisotropic Heisenberg model 

N 

H = J^ [SfSf +1 + S?S] +1 + AS|S? +1 ] (1) 

in the gapped antiferromagnetic region A > 1. For infinite A, excitations are given 
by localized domain walls [6] which can be seen as solitons of unit length in the 
lattice spacing. States involving pairs of such solitons contribute to neutron scattering 
amplitudes 7 J in the vicinity of this anisotropy limit. In the region 1 < A < oo, these 
solitons acquire a finite extent, and their structure and behaviour requires methods 
such as integrability to be properly understood. Experimental realizations of spin 
chains in this regime include CsCoCbj [HI [2] and CsCoBr3 [TOl [TTj . An interesting fact 
is that there exists a clear separation of energy scales between states with 2,4 and 
higher spinon numbers, which provides a way to isolate these separate subclasses of 
states in measurable response functions. We will in what follows be interested in the 
transverse dynamical structure factor (TDSF), defined as 

S- + (k,u;) = - J2 e- ik{j - f) / dte^(Sr(t)S+(0)), (2) 

where the angular brackets (...) denote a zero-temperature average (ground-state 
expectation value). In zero field, which is the case that we will be concerned with, 
we have S~ + (k,w) = S + -(k,uj) = 2S xx {k,u) = 2S yy (k,u). Our purpose will be to 
compute the contribution to this quantity coming from all two-spinon intermediate 
states, using two separate methods relying on integrability. The first method applies 
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to finite lattices, and makes use of determinant expressions for spin operator form 
factors derived within the Algebraic Bethe Ansatz [TH H3] and used to compute 
structure factors of Heisenberg chains for general fields and anisotropies, both for 
two-particle [21 H51 HE] and general multiparticle states [T71 [TS] . The second method 
starts from an algebraic analysis of the infinite chain in zero field [19], and uses the 
quantum group symmetry of the model to express states and form factors directly in 
the thermodynamic limit. This method was first used to obtain the two-spinon part 
of the structure factor for the isotropic (XXX) antiferromagnet in zero field [20l [21] , 
and can also give the four-spinon part [211 [53] . The two-spinon part of the A > 1 
regime was dealt with extensively in an earlier paper |24j . which we would here like 
to revisit with a number of observations. 

The paper is organized as follows. We treat the finite and infinite chains in 
parallel, beginning with a characterization of eigenstates, emphasizing on the two- 
spinon states which will be used to perform the partial trace over intermediate states 
when calculating the structure factor. Section [3J is then devoted to the step-by-step 
construction of the TDSF, gluing together the necessary dynamical constraints and 
form factors. Our results are then presented and discussed, making the correspondence 
between the finite and infinite chains explicit. Sum rules and anisotropy effects 
are then discussed, and we end by offering some perspectives on possible future 
developments. 



2. Eigenstates 



2.1. Finite lattice Bethe Ansatz 

In the gapped antiferromagnetic sector A > 1, the eigenstates of ([T]) are Bethe 
wavefunctions which are each individually characterized by a set of rapidities {Xj} 
solving the Bethe equations 



M 

N ^ 

fc=i 



3 (A J -A fe ) = 2 7 r^. 



(3) 



Here, N is the number of sites (which we take to be even) and M is the number of 
reversed spins starting from the ferromagnetic reference state with all spins up (the 
magnetization is thus S* ot = N/2 — M). Ij are quantum numbers (half-odd integers 
for M even, integers for M odd) which uniquely specify a set of rapidities (and thus an 
eigenstatc) through the Bethe equations. The kernels appearing in the Bethe equations 
are given by 



j(A) = 2 atan 



tan(A) 



2tt 



A 



(4) 



tanh(n?7/2) 

The integer part of the kernel (second term) guarantees monotonicity for real A. This 
is important for the classification of states, which we will discuss below. 

The energy of a state in an external magnetic field h is given as a function of the 
rapidities by 



e =jj: 



sinh 



V 



coshry 



cos 2A, 



4 



M), 



(5) 
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whereas the momentum, which is the sum over quasi-momenta associated to each 
rapidity, has a simple representation in terms of the quantum numbers: 

E' 1 1 TsinfA, + i??/2)l 2tt , „ 

- In -r±r = ^ M m od 2tt. 6 



N 

3 = 



The ground state is given by J? = — + j, j = 1, ...,M, and is not degenerate. 
However, there exists a quasi-degenerate state with momentum tt, given by adding an 
Umklapp to the true ground state. This state becomes exactly degenerate with the 
true ground state in the thermodynamic limit, whereas all other excited states remain 
gapped by a finite value. 

Excited states are constructed by modifying this choice of quantum numbers. In 
everything that follows, we concentrate on solutions to the Bethe equations ([3]) taking 
the form of sets of real rapidities only, {Xj} with Xj € K Vj. We thus do not consider 
string states; however, the important two-spinon states fall into the class of states we 
consider. 

Due to the monotonicity of the kernels 9 n , the left-hand side of © is also 
monotonic in Xj . We therefore have Xj < Xk iff Ij < Ik • Equal rapidities do not yield 
proper Bethe wavefunctions, and we therefore only need to consider sets of distinct 
quantum numbers. Moreover, in view of the parametrization of momenta in terms of 
rapidities used in ([6|), we can restrict to solutions of the Bethe equations in terms of 
ordered sets of rapidities within an interval of width tt, i.e. such that Am < Ai + tt. 
Substituting this condition in the difference of the Bethe equations for Am and Ai 
gives a constraint on Im — I\ , namely 



2tt , s 2tt 

^(Im-/ 1 ) = 2 7 t--X: 



Am — Ai 1 

7T 2 



A; -A! 1 

TT 2 



<2tt(1- — ) (7) 



so Im — h < N — M iff Am — Ai < it. Denote as Cm the set of all ordered sets 
of quantum numbers such that Im — I\ < N — M. Since A is identified with A + n 
as far as the wavefunctions are concerned, this set does not lead to a one-to-one 
mapping of quantum numbers with wavefunctions. We can define a transformation 
S : {(Xjjj)} {(Xj,Ij)} with (Xjjj) = (X j+1 ,I L+1 + 1) for j = 1,...,M- 1 
and (X m Jm) = {X% + tt, h + N - M + 1). Then, Xj < X k iff Ij < I k Vj, k and 
Im — Ii < N — M. S and its inverse are in fact the only two transformations satisfying 
these properties. 

Consider now the set Ai = {{Ij} \ \Ij — i\ < N ^ > 11 }, i-c. the set of 
width N — M of possible quantum numbers, right-shifted by i. For clarity, A$ 
represents the usual set of quantum numbers used in the XXX or gapless XXZ 
cases to get finite, real rapidities. Consider now an element x £ Ai U Ai+\. Since 
Ii — I\ > 2 and Im > N ^^ 1 ^ 1 + i under the transformation S defined above, 
we have S(x) ^ A t U A i+1 Vx. Taking y e (A t U A l+1 )/(A l+2 U A i+3 ), we have 
S(y) G Ai + 2 U Ai + 3 \fy. Therefore, as far as wavefunctions are concerned, we have 
the equivalences Ai U A i+ i ~ A i+2 U A i+3 . We can thus obtain a single enumeration 
of all wavefunctions with purely real rapidities by restricting to Aq and A%, paying 
attention to double counting. The number of different wavefunctions is thus given 
by the number of elements in Ao and A\, minus the number in the overlap Aq n A\, 
namely 

N-M\_(N-M-1\_ N ( N-M \ 

M I { M I' N-M\ M ■ (8 - ) 
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Figure 1. The set of two- 
spinon excitations over the 
ground state of the A = 
4 chain, for N = 200 
sites. Each point repre- 
sents a state in the set 
of excitations defined as 
Ao in the text. The set 
A\ would lead to a very 
similar picture, shifted by 
7r in momentum. The 
darker points correspond 
to states in the set Aq n 
A\, for which we must 
explicitly prevent double- 
counting when summing 
for the structure factor. 

For the transverse structure factor at zero field, we are interested in states in the 
M = N/2 — 1 subsector. At finite size, the true ground state of the system has zero 
momentum, and the N(N + 2)/8 states in Aq are interpreted as the two-spinon states 
constructed above this ground state (see Figure [J). We can interpret the N(N + 2)/8 
states in A\ as the set of two-spinon states built on the 7r momentum quasi-degenerate 
ground state. These would form a figure very similar to that obtained for Aq, but 
shifted by 7r. The overlap of the A a and A\ sets is simply given by the set of states 
living on the darker line in Figure [TJ Counting states only once over this overlap 
region, the total number of states with real rapidities is thus in this case ^j- for all 
momenta considered together, and ^ for each individual allowed momentum. These 
are the two-spinon states of the finite chain, from which we will obtain the finite lattice 
structure factor later on. 

2.2. The infinite chain: algebraic analysis 

We now turn to a complementary approach, which is valid strictly and only for the 
infinite chain in zero magnetic field. Known as the algebraic analysis, it provides 
a construction of states and form factors from purely algebraic considerations, by 
exploiting the quantum affine symmetry of the model. It is described in detail in [19J, 
which contains essentially all the results we need. 

As compared to the finite chain, the infinite XXZ antiferromagnet in the massive 
regime has a number of important characteristics. First of all, the ground state 
becomes exactly degenerate with the 7r-momentum Umklapped ground state. Within 
the algebraic analysis, this translates into the existence of two vacua |vac)(j) whose 
label i = 0, 1 differentiates two different possible asymptotic conditions, specifying the 
z-component of the spin on a chosen reference site. Secondly, a separable subspace of 
the Hilbert space can be identified, which is spanned by multispinon states [5]. The 
lore is that all physical properties of the infinite chain can be derived using only the 
subspace T of (up to denumerably infinite) spinon excitations over the two vacua. 
This is in fact a remarkable statement, considering that the true Hilbert space is 
only isomorphic to T in the subspace of denumerably finite spinon numbers. This 
subspace is spanned by states |£ m , . . . , £±) e ei .,^ with j — l,...,m spinons, each 
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characterised by a pair with spectral parameter £j G {C : = 1} living on 

the complex unit circle and index ej G {— , +} giving the spinon's spin orientation. 
Spinons are always created in (multiple) pairs by local spin flips, and two-spinon states 
are always four-fold degenerate because of the different spinon spin orientation choices. 

If we let Hoc denote the Hamiltonian of the infinite chain |] and T denote the 
translation operator by one site, we have that 



with 



Hoo |£ m , • • • 



,...,e i; (i) 



E m ({ti})\L 



,£1} 



,6 i; (i) 



(9) 
(10) 



(11) 



e(£j) and p(£j) are respectively the energy and momentum of one spinon, with values 



1 - q 2 d 



where 



6 ? (w) = (w;q) 00 (qw x ; g)oo(g; g) c 

OO 



(12) 

(13) 
(14) 



71=0 



Here, q is the so-called deformation parameter of the quantum group U q {sl2), related 
to the anisotropic parameter by the formula A = — Aoo = — (q + q~ 1 )/2 (—1 < q < 0). 
Note that the translation operator T maps the m— spinon states from one vacuum to 
the other: it is therefore more convenient to work with translational invariant spinon 
states 

1 



■ ■ ■ ,£i;p) e 



V2 



16, 



,6 



,ei;(o) 



(15) 



with p = 0, 7r. In particular, the two translational invariant vacua, i.e. no spinons, 
read 

1 



10) 



V2 
1 

71 



vac 



vac 



(o) 



(0) 



vac 



vac 



(i) 



(i) 



(16) 
(17) 



and correspond to the infinite size limit of the two quasi-degenerate ground states 
of the finite lattice. Translational invariant states are eigenstates of the translation 
operator T, 



= e i(Pm({t})+P 



|6> 



>£i;p) 6 



(18) 



J Unfortunately, there exist different conventions for the Hamiltonian in the literature. In the 
algebraic approach, the usual choice is to use a ferromagnetic exchanged Hamiltonian, e.g. Hao = 
— ■j X]j( f7 J' J J+i ~^ a j a j+i +AooTj"j +1 ), which we follow for convenience when using this approach. 
To make contact with the more typical physically relevant antiferromagnetic form JTJ, one should 
simply put Aoo = —A, and rotate spins by n around the z axis on alternate sites, which shifts the 
total momentum by it in the resulting structure factor. 
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Finally, a resolution of the identity within subspace T can be written in terms of 
translational invariant spinons as 

I= E E ^ E tU^\^,...,^M em ,...^e u ...,e m (^...,UM (19) 
m>Oei,...,£m ' p=0,ir i=l ^' 

where the contour integral is such that the square of the spectral parameter £ 2 covers 
the unit circle once. Note that the measure in the resolution of the identity is different 
from the one given in [TU] and normally used in the literature (see for instance |24[[2T| V 
The reason for this is to keep the momentum of the spinon within its physical range 
of width 7r 5] and to obtain the usual spinon dispersion relation (see for example [4]), 
in accordance with the physical picture given by Bethe Ansatz on the finite lattice. 
The different resolution of the identity stems from the fact that, within the algebraic 
approach, states with spectral parameters £ and — £ should be identified (see equation 
(A.12) in [IS]). 

It is mathematically more convenient to work with the following parametrisation 
for the spectral parameter £: 

£ = ie i7r/3/2X , —K < f3 < K, (20) 

with (3 defined modulo AK, with K = K(k) the complete elliptic integral of the first 
kind with elliptic modulus k (not to be confused with the external momentum, when 
we deal with the finite lattice) . The deformation parameter q will be minus the elliptic 
nome q e 

-q = q e = e~ 1L %- (21) 

with K' = K(k') and k! = y/l — k 2 the complementary elliptic modulus. The spinon 
energy and momentum then read 

7T IK /ttK'\ 

e{[3) = 7dn(/3), p(f3) = am(/J) + J , I^^-smhi^— j (22) 

where dn(a;) and am(a;) are the Jacobi elliptic functions with elliptic modulus k. The 
spinon dispersion relation is then explicitly given by 

ei{p) = 7yi-fc 2 cos 2 (p), 0<p<vr, (23) 

where the spinon momentum p take values in the domain p G [0,7r]. This dispersion 
relation is illustrated by the darker line in Figure [T] (up to a shift of Ik', since the 
second spinon is then sitting on zero momentum), and becomes the well-known sine 
curve in the isotropic limit. 

As it is pointed out in [19j . these formulas correspond to the ones appearing 
for example in [5S1 [55] , whereas one might have thought that the range of the spinon 
momentum was 2-7r in the former approach. Although mathematically equivalent when 
the states are correctly identified and the dynamical constraints correctly solved, we 
prefer to work with the spinon momentum p{[3) confined within the interval [0, it] 
to make contact with the finite N Bethe Ansatz clearer. This explains our choice 
—K < (3 < K and the change of measure in the resolution of the identity (TIT))) . 

In the thermodynamic limit, the intermediate states we will use will be 
constructed by creating a pair of spinons over either of the translationally invariant 
vacua. On the finite lattice, these states then correspond to the two-spinon states 
created on the two quasi-degenerate ground states. Formally, there then exists a 
double counting of states on a line (darker line in Figure [T]) in the algebraic approach, 
but this set of relative measure zero does not influence the result for the TDSF. 
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Let us now turn to the construction of the contribution to the TDSF coming from 
the two-spinon intermediate states. As before, we treat the finite lattice and infinite 
one in turn. 

3. The transverse dynamical structure factor 

The gapped Heisenberg chain, as detailed above, has two quasi-degenerate ground 
states at finite N which become exactly degenerate in the thermodynamic limit 
N — ► oo. The chain thus has a tendency to spontaneously develop a staggered 
magnetization along the z axis, selecting for example one of the states \vac)u-\ with 
broken translational symmetry. We here wish to begin by discussing the effects of this 
tendency to spontaneously order on the TDSF. 

When spontaneous order develops, due to the broken translational invariance of 
the ground state, the period of the magnetic structure is doubled, and the magnetic 
Brillouin zone is therefore halved. States and excitations can thus be written as 
modes over this reduced Brillouin zone (since the Brillouin zone is halved, the number 
of modes per momentum is doubled, i.e. we could separately classify sites into odd 
and even ones). One might expect the TDSF to share this k — > k + ir periodicity, 
but this is not the case. Namely, consider calculating a structure factor on either of 
the two possible ordered states |0) or |1) (which can be here taken to represent the 
finite-size equivalents of |vac)(j)): 

Sffaw) = ifYe-toW dte^iilSmSimi). (24) 



i,r J -°° 

We can separate this into the following contributions: 

Sf(k,u) = Ff' l (k,Lj)+e- lk F^ l {k,u J ) (25) 

with F®,£(k,u>) — F^i^ Q {k,u)) + i* 1 ",^ (fc, cj) being the contributions from correlators 
involving pairs of sites at an even/odd distance, which are further separated in the 
functions 

/oo 
dte^(i\SS n+e (t)S b e (0)\i), 
-oo 

/oo 
dte^(^ n+1+e (t)S b A0m, (26) 
-°° 

with e — 0, 1 labelling the sublattice of the base site. All functions F ab '' 1 are then 
manifestly 7r-periodic in momentum, FV* e (fc+7r, a;) = F^ b £ e (k, w), but Sf b (k,uj) does 

not share that property since the e~ lk term in (j2l)|) does not vanish. The points k and 
k + 7r therefore remain inequivalent in the structure factor. The conclusion remains 
the same for the combination S ab = ±(S$ b + Sf) which we consider at infinite size 
(which thus corresponds to the TDSF at a temperature much higher than the splitting 
of the two ground states (which is zero), but much lower than the gap); in fact, the 
structure factor is the same on either ordered state, so all linear combinations give 
the same result (at finite size, we calculate the true zero-temperature TDSF, namely 
on the true ground state only). The limits N — > oo and T — + commute, and 
spontaneous symmetry breaking does not affect the periodicity in momentum of the 
structure factor, which remains 2n. This is similar to what happens in the presence 
of a staggered field, see e.g. the discussion in [57]. With this in mind, let us now 
explicitly construct the TDSF, first on the lattice, then in the continuum. 
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Figure 2. Square form 
factors coming from two 
spinon excitations over 
the true zero momentum 
ground state (— ) and 
over the 7r-momentum 
quasi-degenerate ground 
state (+), for momen- 
tum k = it/ 4 (lower 
curves) and k = 3n/4 
(upper curves), for A = 
4. 



3.1. Finite chain 

On a finite lattice, the transverse dynamical structure factor ([2|) is defined for discrete 
values of momenta k = 2irn/N with neZ. By periodicity over the full length of the 
chain, we can restrict the momentum to the first Brillouin zone n g {0, N — 1}. Using 
the Fourier transform of the spin operators S% = X)j e_ an d introducing 
a formal sum over intermediate states |a), the TDSF can be written in a Lehmann 
representation as 

S-+(k,uj) = 2^|(0|Sfc \a)\ 2 8{uj~ E a +E rj ) (27) 

a 

where Eq is the energy of the ground state |0), and E a is the energy of state \a). The 
form factors are then computed using the known determinant representations [121 [T5] 
once the Bethe equations have provided the rapidities of both eigenstates involved (see 
e.g. (TTl [TB] for a more extensive discussion). 

The full set of intermediate states {\ct}} for a finite lattice separates into many 
different classes of states, only some of which are important for the TDSF. In the case 
at hand, the full Lehmann sum in (|27p can be truncated very efficiently by considering 
only states with real solutions to the Bethe equations, which correspond to the two- 
spinon states described earlier. The structure factor can thus be understood as sums 
of contributions coming from spinons built on the zero-momentum true ground state 
(the set Aq of excitations), and of those built on the it- momentum quasi-degenerate 
ground state (the A\ set). These two contributions turn out to be equal to one another 
(up to 1/iV corrections). This is shown in Figure [21 and illustrates two important 
things. First, both base states (0 and it momentum ground states) in fact yield 
the same structure factor in the continuum limit, so which linear combination of 
these states is taken is immaterial in the definition of the structure factor. Second, 
the structure factor will clearly be asymmetric with respect to the it/2 momentum 
line. In the Ising limit, the form factors do not depend anymore on energy, and the 
structure factor thus precisely follows the density of states. Since the two sets of form 
factors from two-spinon states are then identical (up to 1/N corrections) modulo a 
shift of 7r in momentum, only then does the TDSF becomes symmetric with respect 
to the k = it/2 line. On the other hand, in the isotropic limit, states in the set 
A\/(Aq n A\) contain one rapidity equal to tt/2 (or oo in the usual parametrization 
at the isotropic point), which represents the action of the total spin lowering operator 
SZ—q at zero momentum. Since the Bethe wavefunctions are highest-weight states of 
the SU(2) total spin rotation symmetry, all the contributions from those states vanish 
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due to the selection rules, leaving only the N(N + 2)/8 states in the set A to provide 
nonzero contributions. The TDSF then recovers its well-known 'Viking helmet' shape, 
vanishing at k — and peaking at k = w, and thus having maximal anisotropy with 
respect to the n/2 momentum line. 

Our results for the lattice TDSF are explicitly plotted together with those in the 
thermodynamic limit in Figure El We postpone discussion of these until later, but 
now look at the calculation of the TDSF of the infinite chain. 



3.2. Infinite chain 

In the continuum, the two-spinon part of the TDSF has already been intensively 
investigated [24] . For clarity and completeness, we here reproduce the outline of the 
derivation, which involves some slight differences with the earlier treatment. We start 
by writing the TDSF as 

S- + (Q,uj) = e ~ lQi / dte^(Sj +1 (t)S+(Q)) (28) 

where (■ • ■) = (1/2) X^p=o w (p\ " ' \p)- To decompose this into the different 

spinon contributions, one simply inserts the resolution of the identity ()19|) 

between the spin operators. Upon defining the form factors = 
(j) (vac| a\ |f m , e m ; • • • ; £i, ei) (i ), the TDSF takes the form 

S-+(Q,w)= S ( -+(Q,c) (29) 



m even >o 



with 



lib. I £71 oil ■ 

J i=l Sl p=0,n 

B^m) =E s w(^»» B U({&) = \ x UttM- aX UuM\ a > (30) 

with o = — e lp . We refer to SZ^(Q,u) as the to— spinon TDSF, or simply TDSF( m ), 
which is simply the contribution to the TDSF coming from m-spinon intermediate 
states. Using the [3 parametrisation, the TDSF( m ) becomes 

_2 / i \ m fK rK 

^ (q. ^ = ^ ( 2^ ) / / • • • d ^ ^ - E ™m)] x 

x ^ ^)[Q+p-P m ({/3})]B (ff) ({/3}). (31) 

p— 0,7T 

Let us now concentrate on two-spinon states. The energy of such states is simply 
given by the sum of the energies of the two spinons, 



e 2 (pi,p 2 ) = Iy/l - Pcos 2 ( Pl ) + Iy/i - Pcos 2 (p 2 ), (32) 

whereas the total momentum of the two-spinon state is simply Q = pi + P2 G [0, 2tt\. 
One can write the two-spinon energy as a function of Q and a parameter A = h (pi — P2) 
with A e [— Min(Q/2, n- Q/2), Min(Q/2, w- Q/2)]. The two-spinon dispersion region 
forms a continuous region 1Z = 1Z- U 1Z+ as in the left panel of Fig. [3l which is 
characterized by a very narrow band near Q = and a broad continuum around 
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Q = 7T. This is the usual two-spinon continuum as it appears in the context of Bethe 
Ansatz (see for instance 0], and recall the discussion in the previous section). Some 
detailed comments are however worthwhile to make. If we define k = (1 — fc')/(l + k') 
and 

w±(Q) = -^-v/l + K 2 ±2 K cos(Q), uo(Q) = ~7~7~ sin(Q), (33) 

1 + K 1 + K 

the lower boundary of region 1Z for Q G [0, 7r] is given by 

f «_(Q), Qe[o,Q«], 

wzo(Q)=< wo(Q), Qe[Q K ,7r/2], (34) 
[ w + (Q), Qe[7r/2,7r] 

where Q K = acos(«). In other words, between Q = and Q — Q K , the lower boundary 
is defined by setting A = so p\ = P2 = Q/2 (here and in what follows, we can always 
interchange p\ and P2, and the solution to the dynamical constraints will always have 
this degeneracy). Between Q = Q K and Q = n/2, it is however obtained by setting 
A = ^acos^. Finally, between Q — tt/2 and tt, it is obtained by setting A = Q/2 so 
Pi = Q,P2 = 0, so we here simply fall back on the spinon dispersion relation (shifted 
by Ik' since the second spinon sits at zero momentum). 
The upper boundary of 1Z for Q £ [0, w] is given by 



- / Ik ' + VI " k 2 cos 2 Q, QG[0,Q c ], . v 

where Q c is obtained from solving the quartic equation 

cos 4 Q c - 4cos 3 g c + 2(4/fc 2 - 1) cos 2 Q c - 4(2/fc 2 - 1) cos Q c + 1 = 0, (36) 
whose solution can easily be found in closed form, 

cosQ c = 1 - 4(2/3)^ 11 *> + (2/3) 2 / 3 U ,> S, (37) 



with S = (y3(32 — 5/c 2 ) — Ofc) 1 / 3 . Thus, the upper boundary is simply defined by 
setting A = Q/2 between Q = and Q = Q c , meaning that pi = Q and P2 = 0. 
Afterwards A becomes 0, so the two spinons share the same momentum. The region 
Q € [w, 2ir] is simply described by taking pi — * ir — pi. 

Within region 7Z + , the two-spinon states are therefore ordered in increasing energy 
with decreasing |A|. In region 1Z-, this ordering is reversed; this region exists due to 
the change of sign of the second momentum derivative of the spinon dispersion relation, 
which is greater than zero for small k but becomes negative at an anisotropy-dependent 
value k c = acos-^==. Convolving two spinon dispersion relations thus produces 
this 'folding' region where 1Z- and TZ+ overlap. This region simply disappears at 
the isotropic point, since the spinon dispersion relation then has strictly negative 
curvature. However, a similar region exists for a generic XXZ chain in a field (see the 
discussion in [2"B]y 

For computational purposes and to make contact with earlier results, it is 
convenient to observe that if the region TZ- is reflected around Q = tt/2, we obtain 
the continuous region depicted in the right panel of Fig. [3l which corresponds to the 
sheet C+ in [24] . Note that this reflection is such that region 1Z- is never overlapping 
with 72.+ , but fits precisely under it. The resulting sheet is such that each point refers 
to a single eigenstate (up to the trivial symmetry p\ <-> p2), in other words solutions 
to dynamical constraint equations are unique within the physical region f3 S [-K, K] 




up to simple permutation pi ^ (3%. This was not the case before, since 1Z- and 1Z + 
had a nonzero overlap. 

The boundaries of the sheet C+ in the interval Q <G [0, 7r] are given by 

Oup(Q)=w_(Q) (38) 

for the upper boundary, and 

10 1 , Q S [tt - Q k ,tt] 

for the lower boundary. As in [24], we will calculate the TDSF using the sheet C+ but, 
as we will explain below, one must be careful when dealing with the energy-momentum 
relations. For later use, we also introduce the sheet C_ as the reflection of the sheet 
C + around Q = tt/2. 

To calculate the two-spinon TDSF, we will need to explicitly solve the two-spinon 
energy and momentum relations 

Q = am(/3 1 ) + am(/3 2 )+ 7 r, j = [dn(/3i) + dn(/3 2 )] . (40) 

We follow the derivation in [!23] using addition formulas of elliptic functions. Upon 
introducing the new variables (3± = (f3i ± /3a)/2, the solution to the set of equations 
([4*0]) in the sheet C + is given by 

/3i +) (Q,c) = -i±^F[arcsin(^), K , 



o(+)/vi \ a -1 I 1 + cos (<2) w 2 -Kwg+T \ 

\ |sm(<y)| y + - r / 



with 



T = T(Q, lu) = - n^Juj? - u 2 . (42) 



Within the interval —K < /3i,/?2 < another solution is naturally obtained by the 
transformation f3\ 02, which keeps the set of equations ([40]) unchanged. 

While we have obtained the solutions within the sheet C+, the two-spinon 
continuum is the region 1Z = TZ— U 1Z+. To find the corresponding solution within TZ, 
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we note from the set of eqs P0|) that solutions for (Q,co) and (n — Q,w) are related 

by 

/?i(Q,w) = -&(*-- Q,w), f3 2 {Q,u) = -fo(Tr-Q,uj)-2K. (43) 

While the minus sign in these expressions is not important when evaluating the B 
functions, the shift in 2K is important, since it changes the relative sign of the spectral 
parameters £1 and £2- Defining 



1 + K 



arcsm 



f3^(Q,uj)= dn" 1 



1 + a cos(Q) / ijj 2 — nut 2 + T 



|sin(Q)| 



klOq — T ' 



(44) 

WQ - X J 

the solutions in the two-spinon continuum 1Z are given by 

((3i,p 2 ) = [(3{ +) (Q,u),f3i + \Q,Gj)] {Q^)eTZ+ (45) 

(l3 1 ,f3 2 ) = {f3 { -\Q,u),f3 { - ) (Q 1 u)-2K} {Q,w)eK-. (46) 

We can now substitute these into the expressions for the form factors to obtain the 
TDSF. Recalling formula (|3"Tj) . the projection of the TDSF onto the two-spinon band 
takes the form 

x 2 rK ,K 



Spj(Q,w) = 



df3 1 d(3 2 S[uj-E 2 ({f3})} x 



K J ~K 



(47) 



2 \2KJ 

x e ^)[Q+p-p 2 (mp w (W)- 

p— 0,7T 

The two-spinon form factors involved in (|47[) appeared in the definition of the weights 
B m equation l|30|) . Due to spin conservation, only the form factor with spin 

orientation (£1,62) = (+7+) gives a non-trivial contribution. This is related to the 
form factor with (-,-) through w (vac^ |£ 2 , £1) — ;(i) =(i-i) (vack =F l6>Ci)++;(l-i)- 
All these form factors are obtained from [TI5], and the fundamental building block 
reads 

(i) (vack+|6,6)-- ;(i) 
= (-9) 1 -^i 1 - i r i (9 2 ;? 4 )oo(? 4 ;? 4 ) 3 00 P 2 

with 

(w;g,p)oo 



7(g/tf) 



© 9 *(-« 2 -V l ) (48) 



of Jacobi elliptic functions, one obtains the following expression for B^\[3±) 

\ dn 2 (/3_) 



so that 7(10) 
of Jacobi elli 

(/3±) -lvJ ^)T 

where we have defined 

^(/3)=exp 



JJ (l-wq n p m ), 

n,m— 

r x _ ((-g) 1+ V^;g 4 ,g 4 )oo((-g) 1+g W - 1 ;g 4 ,g 4 )oo f , 

7<H j - ((-g) 3 +V«;;g 4 ,3 4 )oo((-g) 3 +^- 1 ;? 4 ,« 4 )oo 1 ' 

7_(u>) and p 2 = 7 + (g~ 2 ). Using standard definitions and properties 



- /c 2 sn 2 (/?_)sn 2 (/3 + ) 



[fc'<5 CT ,_+dn z (/3 + )^+J 



^ e ke cosh(2fce) cos(2(3ke/ K') - 1 
"IT 



L fc=i 



sinh(2fce) cosh(fce) 



(50) 



(51) 
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with e = and where $d(f3-) is the Neville theta function. Note that B^ a \f3±) is 
invariant under the transformation f3\ «-> fa- 

From (|47[) . we can see that the TDSF( 2 ) consists of the sum of the two-spinon 
region 1Z weighted by ({/?}) plus the same region shifted by 7r and weighted by 
£?( + )({/3}). By noticing that a shift in 2K in one of the parameters (3 implies a change 
in the weights B^ <-> B^ + \ the preceding description is equivalent of weighting 
the sheets C CT with the weights for er = ± and using as the solution to the 

energy-momentum equations the expressions (|44|) without a shift of 2K. Evaluation 
of B (a \f3±\Q , u)) = Bi a \Q,uj) using such expressions gives 



Upon defining 



2K(k)Y l + acos(Q) i%[/?i g) (Q,uQ] 
^o 2 ^D9^(Q,w)] 



1 — K , 



'2K 

MQ,u) = [ — 

7T 



1 + K L 

W1W2 



dE dP 
W2W1 



2K(k)\ ujTW a 



with ujq = ujq(Q), T 



T{Q,u) and W c 
W a {Q,u3)i 



W„{Q,uj) with 



,008(9) 



and considering the multiplicity of the solutions, we can finally write 



SZ+(Q,u>) = 



J a (Q,Uj) 



(52) 



(53) 



(54) 



(55) 



or more explicitly 



1 1 

2~Jf 



1 



E ■ 

*e{-.+} 



acos(Q) i> 2 A (^ ] ) 



ltd — KUJn, 



(56) 



with I(Q ;LJ )gc„ being one if (Q,w) is within the region C CT . Note that this result does 
not agree with the one in [24] . since there is only one weight per sheet and not the 
sum of both weights. Since the two weights are different, the TDSF is asymmetric 
around Q = n/2 for any value of anisotropy A < 00 away from the pure Ising limit. 

We can also see how in this case the result in the isotropic point naturally arises. 
Indeed, as we can see from eq. (|56|) the isotropic limit A — > 1 corresponds to K — ► 1, 
which implies that the weight in the sheet C+ vanishes and only the weighted sheet 
C_ remains. Conservely, in the Ising limit k — ► 0, this asymmetry disappears. This 
is therefore in complete correspondence with what we have described earlier for the 
finite lattice. 
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3.3. Results 

Let us now present our results for the two-spinon part of the zero-field TDSF in 
the gapped antiferromagnetic regime. In Fig. HJ we plot the TDSF over all values 
of momentum and energy covered by the two-spinon states, for four values of the 
anisotropy parameter. As can be seen, by varying A we smoothly go from the Ising- 
like limit to the one of the isotropic point. The finite size results are not plotted 
here, since they are essentially identical to the ones presented. Figure [5] provides 



A=16 




2.5 10 



1.5 3 




_\=2 




3.5 
3 

2.5 
2 

1.5 
1 

0.5 



Figure 4. Two-spinon transverse dynamical structure factor at zero field for 
A = 16, A = 8, A = 4 and A = 2. The approach towards the isotropic limit is 
clearly seen, featuring vanishing of the gap and of the structure factor at the zone 
edges. The asymmetry around the 7r/2 line is clearly seen to disappear only in 
the Ising limit A — > oo. 



a set of fixed-momentum cuts for the same four values of anisotropy, this time 
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showing both the infinite size (solid line) and finite size (colored points, computed 
for a lattice oi N = 1600 sites) results. The accurate agreement between the two 
approaches demonstrates that we have put the building blocks together in the right 
way. As explained in [Mj, the TDSF is characterized by square- root cusps at the 
lower and upper thresholds of the two-spinon continuum, except for Q within the 
range Q K ,ir — Q K , where the structure factor obtains a square root divergence at 
the lower threshold (the latter being given by ojo(Q)). At the isotropic point, the 
divergence at the lower threshold covers the whole momentum interval. Note that our 
results are however different from those in [53], for the reasons explained above. 

The TDSF obeys a number of sum rules, two of which are of particular importance 
for our purposes. First of all, the total integrated intensity is such that 



Perhaps more importantly, the first frequency moment at fixed momentum obeys the 



Since these appear in the Hamiltonian, the value of the right-hand side of (|58|) is easily 
computed from the ground-state energy and its anisotropy dependence. 

In the Ising limit, two-spinon states completely saturate both sum rules. In 
the isotropic limit, it is known [21] that two-spinon states carry 72.89% of the total 
integrated intensity and 71.30% of the first frequency moment sum rule. Figure [6] 
presents our results for the contribution of two-spinon states to both of these sum 
rules. The quantity 02(A) is defined as the fraction of the total integrated intensity 
carried by two-spinon states, and similarly (72 (Q, A) is the fraction of the first frequency 
moment at fixed momentum carried by two-spinon states. 

The total integrated intensity quickly becomes saturated to high accuracy when 
A goes deeper in the gapped regime. This is plotted in the first panel of Figure [6] 
A more interesting point is that (72 (Q, A = 1) does not depend on momentum |21j, 
but develops such a dependence away from the isotropic point. This is illustrated in 
the second panel of the Figure [6] Most of the momentum dependence occurs close 
to the Q = 7T point: in the isotropic limit, the right-hand side of |58|) vanishes for 
Q = it (and so does the TDSF), but not for A > 1. As a function of A, an interesting 
non-monotonic structure is seen as a function of momentum, which is illustrated in 
the third panel of Figure [S] The two-spinon intermediate states clearly carry the bulk 
of the TDSF in this regime. 

4. Conclusions and perspectives 

The understanding of the dynamics of strongly-correlated systems is clearly one of the 
most challenging and long-standing problems in condensed matter, and integrable 
models now provide a pathway towards achieving this goal. Interestingly, two 
independent treatments can be offered in the zero field chain case we considered, 

§ Remember that we have used different definitions of the Hamiltonian for the finite and infinite 
lattices. In the structure factor expressions, k on the finite lattice is thus equivalent to tt — Q on the 
infinite one. 




(57) 
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A=16 



A=16 




Figure 5. Fixed momentum cuts of the two-spinon transverse dynamical 
structure factor at zero field, for the same anisotropy values as in Figure [4] The 
continuous lines are obtained from the algebraic analysis expressions, and the 
points are obtained by a smoothing of the finite lattice result for TV = 1600 sites. 
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Q 



Figure 6. Left (inset): fraction of the total integrated intensity sum rule carried 
by two-spinon states, as a function of anisotropy. Left, main figure: first moment 
sum rule fraction carried by two-spinon states, as a function of momentum, for 
different values of anisotropy. Right: anisotropy dependence of the first moment 
sum rule fraction for fixed values of momentum. 



based either on integrability of the finite lattice, or on the quantum group symmetry 
of the infinite chain. Both approaches have their advantages and disadvantages: the 
former is applicable to generic chains at generic magnetic fields, but is restricted to 
finite chains. The latter is only valid for zero field, and cannot be used to understand 
finite size effects. Together, however, the two approaches paint a rather complete 
picture for zero field. 

In future publications, we will consider the longitudinal structure factor S zz (k, uj) 
in the same regime, as well as the four-spinon contribution to the transverse structure 
factor, thereby generalizing the recent results of the isotropic case [23] to this sector. 
An interesting further line of investigation would be to apply the algebraic approach 
within the gapless anisotropic regime. Finite temperature results would also be of 
great interest. 
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